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^ • We present the new expressions for the formfactors of local operators for the XX 

- quantum spin chain as a Cauchy determinants. Using the known functional form of 
the correlator at large distances we propose the new expression for the constant for the 
asymptotics of the correlator as a Cauchy determinant. We calculate the momentum 
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, distribution for the general case of the XXZ - spin chain and point out that it is com- 

I pletely different from that for the Luttinger model (the system of fermions). For the XX 

l/^ I chain we compare numerically the value of the lowest formfactor and the expectation 

' value of momentum- zero operators which is determined by the functional form of the 

correlator. 
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^ I 1. Introduction. 

Ph| In the present paper we study the spin-spin equal-time correlator in the XX - spin chain 

I with the Hamiltonian H = ^ Z^i^i (Ci^f^i+i + o'l^'o'f+i), which is the particular case of the XXZ 
J3 I spin chain with the anisotropy parameter A = 0. It is worth mentioning that although the 
XX-chain can be solved with the help of mapping to the free-fermion system via the Jordan- 
Wigner transformation the system is essentially the the hard -core bosons on the lattice (at 
^ ' the half-filling) and is not equivalent to the free fermions which manifests for example in the 
' completely different correlators in the two models. The interest in the calculation of correlators 
for the XX spin chain is in particular due to the possibility to use it as a testing ground for 
various approaches for the calculation of the correlation functions for the general case of the 
XXZ quantum spin chain and the other models solvable with the help of the algebraic Bethe 
ansatz method. The correlators for the XX spin chain as well as the correlators for the Ising 
- type spin chains both at zero and finite temperature where previously studied long time ago 
[1], [2], [3], [4] using various methods based on the application of theorems for the Toeplitz 
determinants. Recently it was also shown [5] that the equal-time spin-spin correlator in the 
XX- model can be calculated exactly on a chain of finite length at any distance x of order of 
the chain length and the functional form of the correlator coincides with that predicted by 
the conformal field theory [6] or the bosonization procedure [7], [8], [9]. 

The goal of the present paper is twofold. First, we study the consequences of the functional 
form of the correlator for the XX chain in particular for estimates of the certain formfactors 
of local operators for this model. The interest in the formfactors is in particular, due to the 
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existence of the closed expressions for the formfactors for the XXZ model in the form of the 
determinants obtained in the framework of the algebraic Bethe ansatz method (for example, 
see [10]). Here we present the new expressions for the formfactors of local operators for the 
XX - quantum spin chain as a Cauchy determinants. Using the known functional form of the 
correlator at large distances we propose the new expression for the constant for the asymptotics 
of the correlator as a Cauchy determinant or in the form of the product depending on the 
momenta. At the same time using the simple scaling arguments the numerical value of the 
"lowest" formfactor is found. This value is compared with the values of the other formfactors, 
which shows their drastic decrease with the increasing of the energy or momentum of the 
intermediate state. We discuss the consequences of this fact. The second goal is to calculate 
the momentum distribution for the general case of the XXZ - spin chain and point out that 
it is completely different from that for the Luttinger model (the system of fermions). The 
spectral properties (momentum distribution) are of interest for both models due to their 
possible realizations in nature. In particular we point out the universal character of the 
connection of the momentum distribution singularity with the asymptotics of the correlators 
and the constants in front of the asymptotics. For the Luttinger model the connection of the 
momentum distribution with the functional form of the correlators is obscure in the literature. 
The momentum distribution for the XXZ spin chain have not been considered previously. 

In Section 1 we show how the usual bosonization procedure leads to the functional form 
predicted by the conformal invariance and discuss the momentum distribution for the system 
of spinless fermions (Luttinger model). In Section 2 we calculate the formfactors and obtain 
the new expression for the formfactors in the form of the Cauchy determinant. In Section 
3 we briefly review the well known calculation of equal-time spin-spin correlator in the XX 
-model [1], [2] and present the expression for the correlator at the distances x ~ L. We then 
use the functional form and the expressions for the formfactors to get the new expression for 
the constant in front of the asymptotics. Later in this section we evaluate the momentum 
distribution for the XXZ spin chain. In the Appendix we present for completeness some 
formulas and the numerical constants used in the text and also briefly review the previous 
calculation of the constant in front of the asymptotics for the XX spin chain. 

1. Bosonization. 

Consider the effective low -energy Hamiltonian which build up from the fermionic operators 
(ofe, Cjfc, k — 2Trn/L, n & Z, L- is the length of the chain) corresponding to the excitations 
around the right and the left Fermi- points and consists of the kinetic energy term and the 
interaction term H — T + V with the coupling constant A: 

H^Yl K(4(^k - 4ck) + \/L a^ak+g4,Ck'-g. (1) 
Deflning the operators [7] 
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where |/c + p| < A, where A is some cut-off energy, which for the states with the filled 
Dirac sea have the following commutational relations 

[Pi(-P); Pi(p')] = ^^P'P' [P2{p)iP2{-p')] = ^^p,P'^ 
one can represent the Hamiltonian in the following form: 

H ^-j-Yl (Plb)Pl(-p) + P2{-P)P2{P)) + A 53 — (Pl(p)p2(-P) + Pl{-P)p2{p)) ■ 
^ p>0 p>0 

To evaluate the correlators in the system of finite length and make the connection with the 
conformal field theory predictions, one can introduce the lattice fields 711^2(2^) corresponding to 
the Fourier transform of the operators pi,2(p) and represent the last Hamiltonian in the sector 
with the total number of particles and the momentum AA^ = AA^i-|- AA2, AQ = AA^i — AA^2, 
where AA^i 2 are the numbers of particles at the two Fermi-points, in the following form (this 
was first proposed in ref.[12]): 

H = 27rE (^K(^) + nl{x)) + Xn,{x)n2{x)) + ^u{X) \i{ANf + (l/0(Ag)^l , (2) 



X 



where the parameters u{\) = (1 — A^)^/^ and ^ = ((1 + A)/(l — A))^/^. Calculation of finite 
- size corrections to the energy of the ground state for the XXZ- spin chain (see for example 
[13]) leads to the expression (2) and allows one to obtain the parameter ^ which leads to the 
predictions of critical indices according to the conformal field theory. The calculation gives the 
value ^ = 2(71 — r))/ 71, where the parameter 77 is connected with the anisotropy parameter of 
the XXZ - chain as A = cos{ri) (in the present paper we consider the antiferromagnetic XXZ 
chain in the massless regime < A < 1). Rescaling the variables as ni^2{x) = v^27r ni^2{x) 
one finds the commutational relations [ni{x); ni{y)] = d'{x — y) in the continuum limit which 
allow one to introduce following conjugated field and the momenta: 

7r(x) = -^{ni{x) - n2{x)); d^(f){x) = -^(ni(x) + n2{x)) 

4>{x) = N{x) = N^{x) + N2{x), 7Vi,2(x) = r dy h,M 

Jo 

In terms of these variables the Hamiltonian density takes the following form: 
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H = -u{X) il/0^\x) + mi^)r = ^«(A) r{x) + (90(x))2 , (3) 
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where 



7r(a;) = ^(x), 0(x) = (l/v/O0(^)- (4) 

The last equation (4) is nothing else but the canonical transformation, which is equivalent to 
the Bogoliubov transformation for the original operators pi,2(p)- Next to establish the expres- 
sions for Fermions one should use the commutational relations [a'^{x); pi{p)] = — e*^^a"'"(a;) 
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and the same for c^(x). Note that these last relations were obtained using the expression 
with original lattice fermions: = e'^^a"'"(y)a(7/). In this way we obtain the following 

expressions for fermionic operators: 

a+(c+)(x) = exp I ±^ P^e-'A = K,,^ exp {^i27^N,,2{x)) , (5) 



where the fields Ni^2{x) differ by the normalization v27r from the fields Ni^2{x) and Ki^2 are 
the Klein factors - the operators which creates the single particle at the right (left) Fermi 
-points (we omit here the usual exponential suppression e""^'^'/^ in the exponent and and the 
constant factor in front of the exponent 1 / v^tto which in the limit a — > leads to the correct 
anticommutational relations [9]). Note that the above expressions (5) are equivalent to the 
known "field-theoretical" bosonization formulas [8]. 

Now let us apply the above formulas to the specific case of the XXZ - spin chain. Using the 
Jordan- Wigner transformation cr^ = a+exp(i7riV(a;)), where stands for the "original" lattice 
fermionic operator, and performing the obvious substitutions N{x) — > x/2 + Ni{x) + N2{x) 
and — > e'^^^a"^(x) -|- e~'^^^c"'"(x), pp — 7r/2, we obtain after the canonical transformation 
(4) the expression for the spin operator which determines the leading term in the asympotics 
of correlator for the XXZ -chain: 



< ~ (-l)-exp ^-tn^i{N,{x) - N2{x)) j , (6) 

where Ni^2{x) - are corresponds to the free fields 7t{x), (f){x), obtained after the transformation 
(4). To these operators correspond the new operators pi,2{p) and the new fermionic operators 
(quasiparticles) . Analogously the term responsible for the subleading asymptotics has the 
form 

exp [i2T,{l/fO{N^{x) + N2{x))^ exp (-z7r^(7Vi(a;) - N2{x)) ^ 

Averaging the product of exponents in bosonic operators for the expression (6) and using 
the properties of pi,2(p), (pi(-p)pi(p)) = ^0{p) and {p2{p)P2{-p)) = |f ^(p), we get for the 
correlation function G{x) = (0| cr jt^^(7j~|0) the following sum in the exponent: 

where C - is some constant. Then using the formula J2'^=i — — ln(l — z) and substituting 
the value ^ = 2(7r — r^)/?! — > 1 we obtain the following expression for the XX - chain: 

G{x)^Co ,^ a^)-^—^^l/2 {x»l). (7) 
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Thus, although bosonization, which deals with the low-energy effective theory, is not able to 
predict the constant before the asymptotics, the critical exponent and the functional form are 
predicted in accordance with conformal field theory. 
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Let us show how the general form of the correlators, predicted by CFT, including the higher 
order terms in the asymptotics and their functional form can be obtained in the framework 
of bosonization for the XX chain and the spinless fermion model. For the XX chain we 
use the same formula for the operator — a'^exp{inN{x)). Next, the lattice fermionic 
operator should be projected on to the effective continuous operators a'^{x), c^{x) which 
enter the effective low-energy theory (Luttinger model). Clearly, each of the operators should 
be accompanied by the combinations of the type {c'^{x)a{x))'^, {a~^{x)c{x))"^ conserving the 
total number of particles. Then one should consider the following substitution for the lattice 
operator 

m 

where pp — 7r/2 and Ci(m), C2(m) are some constants for an arbitrary integer number m. 
Clearly, repeating the procedure leading to the expression (7) we will obtain the expression for 
the correlator as a sum of the terms of the form (7) with the critical exponents a{m) = a + 
m^/ a in agreement with predictions of the conformal field theory. Note that we considered only 
the operators which conserve the number of particles, in general we would obtain the spectrum 
of primary operators labeled by two quantum numbers n and m of the form a{n, m) = r?a + 
vr^ ja in agreement with the prediction of the Gaussian model of the conformal field theory 
(\/^ plays the role of the compactification radius). In general there can be the contributions 
corresponding to the descendant operators, which amounts to the presence of the additional 
terms of the type (d^)'^ Ni^2{x) in the last equation for aj, however, one can argue that these 
terms are absent in the case of the XX spin chain. Let us stress that the correlator G{x) is 
represented as a sum of the terms, each of the functional form (7) with the corresponding 
critical exponents Q;(m) = a + vn? ja. Note that in the similar way the asymptotics and the 
functional form of each term in the expansion of the density-density correlation function can 
be found. For this case the coefficient for the leading term is equal to unity. 

In Section 3 we will compute the momentum distribution (nfe) = (o"^a^), which is con- 
nected with the asymptotics and the functional form of the correlator G{x) and the behaviour 
of the formfactors at small energies or momenta q. Here before studying the momentum dis- 
tribution and their relation to the formfactors in the XX spin chain, let us comment on the 
momentum distribution for the system of fermions, or for the Luttinger model (1). As an 
underlying model one can consider, for example, the exactly solvable system of fermions with 
the XXZ - type interaction, H = J2i{^\{cit (^i+i + h.c.) + Ariiriij^i) , for which the parameters of 
the Luttinger model and the critical exponents can also be calculated exactly. For this model 
the bosonization procedure leads to the following asymptotic expression for the correlator 
Gf{x) = {a^tto) analogous to the correlator (7): 



sm{ppx) 
(Lsin(f))' 



GH^)=^(PF,A)^^^^^^, (x»l), (8) 
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where pp is an (arbitrary) Fermi momenta, C is some constant depending on pp and the other 
parameters of the model (A) and the critical exponent equals 

a{X) = 1(^ + 1]= ^ > 1. 

Clearly at A = we have A = 0, q;(A) = 1 and the correlator (8) reduces to the free-fermion 
correlator. One can see from (8) that the form of the correlator and the critical exponents in 
the two models are completely different although the thermodynamic quantities are the same. 
Calculating the Fourier transform of the correlator (8) taking into account its functional form, 
one readily obtains the well known momentum distribution for the Luttinger model close to 
the Fermi points: 

(rik) = (a+a,) = ^ + C\k-pp\\ 7 = «(A) - 1 = (^^/^ - (9) 

First, both this expression and the critical exponent 7 are exact in a sense that they are 
not an artifact of approximation of the initial lattice spinless fermion model by the Luttinger 
model. Both the function (9) and the value of 7 are the direct consequence of the equation 
(8) valid in the initial model so that 7 is directly connected with a (A) in eq.(8). Second, 
the constant C in the equation (9) related in a simple way with the constant C{pf,A) in 
eq.(8), the fact that have not been pointed out previously. In fact, the subleading terms in 
the correlator lead to the contributions to (9) which are of higher order in \k —pp] and can be 
neglected in the vicinity of the Fermi point. Since eq.(9) was obtained from (8) for k — > pp, 
the constant C is proportional to C(pf, A) with the coefficient which depends only on the 
exponent a (A) in exactly the same way as for the constant in front of the asymptotics of 
G{x) for the XXZ spin chain. Wc will derive this relation below in Section 3. Although the 
behaviour of the momentum distribution for the XXZ spin chain which is determined by the 
asymptotic behaviour (7) is completely different from that in the Luttinger model, the leading 
order singularity is also unambiguously predicted by the exponent a and the constant before 
the singularity is related in a simple way to the constant Cq. 

2. Formfactors for the XX -spin chain. 

Let us evaluate the expressions for the formfactors of the spin operators af" for the XX- 
model at zero magnetic field with the Hamiltonian: 

^ i=i 

where the periodic boundary conditions are implied. Performing the well-known Jordan- 
Wigner transformation: 

<7+ = e-^(-)a+ = exp(z7r$:nOa+, 

Kx 
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we obtain the following Hamiltonian written down in terms of the Fermionic operators: 



where M is the number of particles. We assume for simplicity L- to be even and M = L/2 
to be odd {S^ — for the ground state and (M — 1)- even, we also assume L to be even so 
that the ground state is not degenerate) . Then the Hamiltonian is diagonalized with the help 
of the following Fourier transform in the different sectors: 

The formf actor of the operator cr^ at the L-th site (which can be denoted also as Uq) takes 
the form: 

(wko iM) = (oi ( -^^""^ (fj^"! ^10) 



where the sets of the momenta {p} = {pi, . . .pu} and {q} = {qi, . . . qM-i} correspond to the 
eigenstates in the sectors with different number of particles (determined by the integers or 
half- integers numbers) . The set {p} corresponds to the ground state of the system while the 
set {q} corresponds to the arbitrary excited state in the sector — —1. One can represent 
the formfactor as a determinant using Wick's theorem and the following expression for the 
average: 

(0|^.<|0) = Z (i_e^(.-.)) = L sin((p-g)/2) - ^^'^ 
Then we obtain the formfactor as a sum of determinants in the form 



^ ^1 sin((p^)-5,)/2), 

where the set of the momenta {p*^™-'} is obtained from the set {p} by the exclusion of the single 
momenta Pm and q = J2iLi^ Qi- Thus we obtain the following expression for the formfactor: 

(Mko~IM) = -^detij{Mijip,q)), i,j = 1,2, . . .M, 

where M x M matrix equals: 

^ T ■ (( ^ MOV J = '^,---M-l, 
L sm{{pi - qj)/2) 

MiM = e-'^''\ i = l,...M. (12) 
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This determinant can be calculated in the following way. Introduce the new set of the momenta 
W} — {^1' • • • Qm-1, Qm} and take the limit qm = iQ, Q ^ oo. Then 



Lsm{{pi-iQ)/2) L 
and the formf actor is represented as a Cauchy determinant: 

detij{Mij{p, q)) = ^e^det^j ' 



2i ''\Lsm{{p,-q'^)/2)J- 

Using the well known formula for this determinant and taking the limit Q — > oo, we get the 
following formula for the formf actor: 

^^^'^^"VlUJ ) n,,sin((p,-,,)/2) ' ^^^^ 

where the relation sin((pj — iQ)/2) —>■ ^e'^*e*^ was used and the products are over the initial 
sets of the momenta {p} and {q}. Let us stress that the expression (13) for the formfactor 
is valid not only for the ground state configuration of the momenta {q} but for an arbitrary 
excited state characterized by the momenta qi, ■ ■ ■ qu-i- The advantage of the expression (13) 
for the formfactors in comparison with the determinant expressions is that using the expression 
as a product one can hope to extract the behaviour of the formfactors on the momenta q for 
the low-energy excited states (for example, the one-particle - one-hole excitations) and thus 
extract the asymptotic behaviour of the correlator. Note, that the expression (13) is valid for 
an arbitrary filling fraction (arbitrary M) or for the XX spin chain in the magnetic field. 

Let us comment on the expressions for the formfactors obtained in ref.[14]. One can use 
the explicit expressions for the wavefunctions for the hard-core bosons in the coordinate space 

• i<j p 

where Xi is the coordinates of spins, e{x) - is the sign function and P - is the permutation of 
M particles, P G Sm, obtained as a limit of the Bethe wave function for the XXZ - chain and 
using the formula 

^e(x)e^(^'-«)^ = ictg((p - q)/2) 

X 

to derive the following formula for the formfactor for the same sets of the momenta {p} and 

({?}ko~l{p}) = deti,(Mi,(p, g)), i, j = 1, 2, . . . M, 
where M x M matrix Mjj equals: 

M,, = ictg (^^) , j = l,...M-l, 
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M,M = 1, ^ = 1,...M. (14) 

Let us show the equivalence of this expression with the expression given by eq.(12). The 
expression (11) can also be represented as 

(0|c,a+|0) = -^(ctg((p-g)/2) + z). 

Then we obtain from (10) the determinant of the matrix of the same form as (14) with the 
matrix elements Mij for j ^ M replaced by — ictg((pi — gj)/2) + 1. To show the equivalence 
of two expressions for the formfactor one can use the following theorem. Consider the deter- 
minant of the sum of two matrices detij(My + Oy), where the second matrix = Cj0j - is 
the matrix of rank 1 . Then we have: 

M 

detij {Mij + ttij) = detij{Mij) + ^ Cfcdeti,(M^^''^), 

fc=i 

where the matrices M^^^ differ from the initial matrix M^j only by the substitution of its k-th 
column by (pf. 



Mi;> = (1 - 5,k)M,j + 



Applying this statement in the case of the matrix aij = 1 (c^ = 1, 0j = 1), to the determinant 
obtained from (10) one readily get (apart from the irrelevant phase factor) the expression 
(14). Thus the equivalence of the formfactor (13) with that obtained in [14] is shown. 

Using the expression for the formfactor (13), it is easy to obtain another expression for the 
formfactor corresponding to the ground-state configuration {p} = {pi, . . .Pm}, Pi — 2t:/ L{i — 
(M + l)/2), (M- is odd), and the ground state in the sector with M — 1 particles {q^} — 
■ ■ ■ ?M-i}, — 27r/L(i— M/2). Equivalently in eq.(12) one can take the shifted momenta 

p, = (27r/L)(z), ^ = 1,...M, = (27r/L)(j + l/2), j = l,...M-l. 

and introduce an extra momentum qu — (27r/L)(M-|-l/2) into the products of eq.(13) in order 
to make the correspondence with the expressions for the building blocks for the correlator Rj^ 
introduced in the next section more transparent. In fact, one can use the simple formulas for 
the product of sinuses [15] 

M-1 M-1 

n sin(^(A;)) = 2(^-^)/2x/Z, [] sin(^(A; + 1/2)) = 2(^-^)/2, 
fc=l ^ k=o ^ 

where M = L/2 (L is even) is taken into account, to obtain instead of (13) the product with 
equal numbers of variables Pi and qj: 

,l,(fnm^- 1 f'\''~' n.<,- sin((p, - p,)/2) U,<j Midi - Qj)m .... 

n,,sin((p.-g,)/2) ' ^''^ 
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where the new sets of the momenta {p} and {q} are: 



= (27r/L) (i) , = (27r/L) (i + 1/2) , i^l,...M. 

In the equivalent form the expression (15) can be again represented as a Cauchy determinant 
of the M X M matrix: 

^»^'°'»-;^^<"( Lin(i(4-l/2)) )- ''-'r (^«* 

This equation for the formfactor will be used in the next section to compare the constant Co 
for the asymptotics of the correlator G{x) with the square of this formfactor (16) and obtain 
the new expression for the constant Co as a Cauchy determinant. 

Let us show that both expressions for the formfactors (12) and (14) leads to the same 
expression for the correlator as an M x M determinant [14], which in turn is equivalent to 
the exact expression for the correlator, first obtained in rcf.[l]. We consider the square of 
the formfator (14) and sum over the intermediate states {q} with the weight e*''^, where 
Q = YjfLi^ li- Since the formfactor is an antisymmetric function of the momenta g,, the sum 
over each q^ can be extended to the whole region q e (— 7r;7r) independently of the other q^. 
We also represent by the single determinant the product of the determinants of two matrices. 
One can use the following formula for the matrix elements (14), 

ctg(^— )ctg(^^) = ctg(^-^) (^ctg(^— ) - ctg(-^) j - 1, 

and the following Fourier transform, 

\ Y: e^^-ctg((g - p)/2) = i(l - 5,,o)e'^^. 

Thus, using the formfactors in the form of the double sum of (M — 1) x (M — 1) determinants 
(with alternating signs and the sets of the momenta {p^^^} which is obtained from the set {p} 
by omitting the momentum p„) one easily obtains the following expression for the equal-time 
correlator as a difference of two M x M determinants depending on the distance x: 

G{x) = detSf ) (M,, + R,j) - detSf ) (M,,) , (17) 

where the matrices M^, Rij are given by the equations: 

L> '''^L tg{{p,-pj)/2) ' ^^:>-L^ 

which was presented in ref.[14]. From the determinant expression for the correlator (17) one 
can obtain the exact expression for the correlator G{x) on a finite lattice which is obtained in 
the next section with the help of the method proposed in ref.[l]. 
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3. Exact correlation function: momentum distribution. 



Let us briefly review the exact calculation [1] of the spin-spin equal-time correlation func- 
tion (density matrix) for the XX- spin chain on finite lattice of the length L [5] : 

G'(x) = (OK,ar|0). 

Using the Jordan- Wigner transformation relating spin operators to the Fermi operators {af , ai 
= exp(i7r n;)a^, the correlation function G{x) can be represented as the following av- 
erage for over the free-fermion ground state: 

G{x) = (O|a+e'"^(^)ao|0), 

where N{x) — Z^fjf'^ni. Introducing the operators, anticommuting at different sites, 



„i7rni 



where rii — aftti - is the fermion occupation number, with the following correlators with 
respect to the free-fermion vacuum, 

(0|SiA,-|0) = 2G'o(i-j), (0|A,A,|0) = 0, (0|5,Ej|0) = 0, 

where the free-fermion Green function on finite chain Gq{x) is 

^ / N ,1-1- , V sin(7rx/2) 
Go(x) = (0|ai.a.|0) = 

one obtains the following expression for the bosonic correlator: 

G{x) = ^(0|5o(^ii^i)(^2i?2) . . . {A,-iB,-Mx\^)- 

Note that we assume the periodic boundary conditions for the initial spin operators, so that 
strictly speaking, the above formulas are valid only for the case when L - is even (the ground 
state in not degenerate) and M = L/2 - is an odd integer (so that L/4 - is not an integer). 
That is easily seen from the boundary term ~ exp(i7r(M — 1)) = +1 when M- is odd, momenta 
of fermions - are integers (not half-integers) and their configuration is symmetric around zero. 
In that case the free- fermionic function Gq{x) - is given exactly by the above formula, (as 
for the case M - even, it is clear from the general grounds that in this case Eq.(9) will be 
modified by the terms of order 1/L for any x). Using Wick's theorem we obtain the following 
determinant oi x x x matrix: 

G{x) = detij{2Go{i - j - 1)), ij = l,...x. 

Due to the form of this matrix {Go{l) — for even I) this determinant can be simplified and 
the following formulas are obtained: 

G{x) = J(i^iv)^ {x = 2N), G{x) = -IrnRn+1, {x = 2N + 1), 



11 



where we denote hy Rn the following determinant of the N x N- matrix: 

Rn = det,,- {i-iy-^2Goi2i - 2j - 1)) , i,j^l,...N, 

where ^0(2;) is the same Green function of free fermions as above. Since Rn - is the Cauchy 
determinant one can obtain the following expression for it on the finite chain: 



TT/ j^J-^ Vsin(7r(2A; + 1) /L)sin(7r(2A; - 1) /L) ^ 



Note that we obtained the exact expression for the correlator on the finite lattice [5]. Prom 
the expression (18) it is easy to obtain the correlator in the thermodynamic limit (L — > 00) 
which is given by the similar product. For finite chain it is easy to evaluate (18) numerically 
and compare the result with the asymptotic (7) at x >> 1 and x ~ L. At the distances 
X ~ L the correction to the asymptotic formula (7) behaves like ~ 1/L. We find that the 
exact correlator coincides with the correlator given by (18) with very high accuracy up to the 
very small distances a; ~ 1. 

Using the expression (18) one can get the new expression for the constant Cq for the leading 
asympotics for the corellator G{x) (7). For convenience one can redefine the spin operators 
according to — > (—1)^(7^, which is equivalent to the change of the sign for the Hamiltonian 
of the XX chain. That leads to the shift of the singularity in momentum distribution from 
the values q = ±n to g = 0, which will be convenient later in this section for evaluating the 
momentum distribution. Then the correlator (7) takes the form: 

G{x) = Co 



(Lsin(7rx/L))°' 



{a = 1/2). According to this formula the correlator G{x) at the distance x — L/2 equals 
G{L/2) — Go/y/Z with the corrections of higher order in 1/L. On the other hand, since our 
expressions for the correlatar are obtained for M — L/2 - odd, we have G{L/2) — ^R m-i R m+i , 
where Rn is given by eq.(18), or, equivalently, 

Rn = detij [ J. — 777^ ) > hj^^,---N. 

\Lsm {[2Tr/L)[i -j- 1/2)) J 

Clearly, this determinant can also be represented as the product of sinuses. First, this relation 
can be used to obtain the new expression for the constant Cq. Second, one can use the relation 
G{L/2) = Cq/s/L to predict the asymptotic behaviour of the last dererminant Rn at large 
L (or M). Next, one can compare Rn with the expression (16) for the formfactor for the 
intermediate state corresponding to the ground state in the sector with M — 1 particles. Note 
that there is an obvious similarity between the expressions for two different quantities. In 
particular using these arguments one can predict the asymptotic behaviour of the formfactor 
■^({9*"'^''}) large L including the constant proportional to Co which is interesting by itself 
since in general it could allow to relate the constant Co to the square of the lowest formfactor. 
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In fact, neglecting the terms of order 1/L, one gets the relation 



where Rm/2 equals to the determinant which is expressed through M as 

^"/' = '"'*-'( Msm((VMK.-,-l/2)) )- *.J = l.-^/2- («) 

Taking this expression into account one can find both the scaling behaviour (the behaviour 
in a power of L) and the constant for the determinant in eq.(16) determining the the square 
of the lowest formfactor In fact, using the scaling arguments, i.e. neglecting the 

corrections of order 1/L, comparing the expression (19) with the determinant (16), one finds 

mW''})f = ^- (20) 

Thus we found the "lowest" formfactor i.e. the single formfactor, corresponding to the ground 
state in the complete set of the intermediate states. This formfactor can be compared with 
the Fourier transform of the correlator G{x) at zero momentum, J2x=i ^(^)) which gives the 
sum of the squares of the formfactors for all intermediate states with zero total momentum 
(including the ground state formfactor (20)). First, one observes that in order to calculate this 
quantity it is sufficient to use the leading order expression (7) for the total correlator G{x). 
In fact, it is easily seen that the sum for the subleading terms is suppressed by the powers of 
L. Second, one can see that to calculate the sum Y.x=i ^(^) sufficient to replace the sum 
by the corresponding integral since the corrections are again of order ~ 1/L. Denoting the 
corresponding sum of the formfactors as I'^oP = (1/-^) ^x=i ^(^) '^^ obtain 

|V,„r = ^, C.^C„l'ay^-±—, (21) 

which clearly has the same scaling behaviour as eq.(20). The constant Ci in the last equation 
should be compared with the constant V2Co in (20). Calculating the integral, one finds the 
coefficient Ci = A/7rr(l/4)/r(3/4)Co = 1.6725. .Co which is to be compared with v^Cq. One 
finds the striking coincidence of two different values, which can be explained only by the 
extreme smallness of the contributions of the excited states formfactors to the sum (21) in 
comparison with the "lowest " formfactor (20). In our opinion this result is of interest in 
connection with the study of the correlators for the general case of the XXZ spin chain. The 
smallness of the higher states contributions in (21) can also be seen from the representation 
(13). For example, consider the configuration {q} with g = with qi and qM-i shifted by one 
step to the left and to the right respectively. One can easily calculate the value of (13) which 
is lowered by the factor 1/4 with respect to the ground state. Taking the square of this value 
one gets the factor 1/16 which is close to the relative difference between (20) and (21). 
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Let us turn to the study of the momentum distribution for the XX spin chain which turns 
out to be completely different from that for the system of fermions (Luttinger model) (9). The 
results obtained are equally well apphcable to the general case of the XXZ chain provided the 
constant for the asymptotics (7) for XXZ case is known [11]. Turning back to the sum of 
the formf actors j'^oP (21) one can see that this quantity is nothing else but the probability 
n{q — 0) to have the total momentum q — which is defined by the equations: 

n(q)^W^~''^G{x), G{x)^^e''^MQ), ri{q) ^ {a+a^), (22) 

^ X q 

where the sum over q = 2TTn/L, n E Z is extended over the interval q G (— vr, vr), and (t~ 
are the Fourier transform of the initial spin operators and for the half-filling the normalization 
condition is I]gn(g) = 1/2. The value of n(0) = Ci/a/L was already calculated above (21) 
(the useful formulas which allow one to perform the numerical estimates are presented in the 
Appendix, note that the value of n(0) can be easily obtained for the XXZ chain). Let us 
stress once more that the value [■^ol^ = ?T'(0) (= Ci/y/L) is an exact value of n{0) - the sum 
of the formf actors with q — in the limit of the large chain L >> 1. Next, from (7) one can 
calculate n{q) for the XXZ chain for the sufficiently small q in the form: 

, , Co ,7ra,2"r(2 - a) r(n + a/2) 27m , , 

Evidently, at g = one obtains the constant C\ calculated above. This formula is correct 
provided the sum over x in eq.(22) can be replaced by the integral i.e. for n « L since at 
larger values of q the oscillatory behavior should strongly suppress n{q) up to the value ~ 1/L 
at ~ ±7r. The behaviour in n can be found using the formula for the ratio of two Euler's 
gamma- functions in the equation (23). For the XX chain it predicts the behaviour of order 
atl<<n<<L where n — Lg/27r. Clearly, at n ~ 1 or g ~ 27r/L this formula 
is not valid, so at n ~ 1 one should use the initial formula (23). In the general case of the 
XXZ chain the equation (23) gives 1/L"n^^". The constant in front of this asymptotics for 
n{q) is not equal to the constant Ci calculated above. However, it can be easily found for the 
XX chain from eq.(23) as C2 = (2/7r)r(3/2)Co = (Iv^)Co = 0.564.. Co. As in the case of the 
Luttinger model, the contribution of the subleading terms in C(x) arc suppressed by a powers 
of 1/L in the region, where the asymptotics 

„(,) = = ^M_L (24) 

is valid i.e. in the region where the value of n{q) is parametrically larger than 1/L. In this 
sense the momentum distribution given by the equation (24) is exact. The constants C2, C'2 
for the XXZ chain equal 

. na 2mi-a) 

C2 a = sm — — —Co, 

z 7r(l — a) 
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and G^ia) = C2(a)(27r)i-". For the XX spin chain we get C2 = 0.564.. Co, C'^ = 1.410.. Co- 
There is no reason to expect that the "singular" part of the function n{q), n{q) = CyLq^^''^ 
(24) saturates the sum rule 

TT ^ 

^(^) = 2' 

since in general the subleading terms give the contribution of order L{l/L) ~ 1. However one 
can calculate the contribution of the function (23), Z^^^o'^l?) taking into account the known 
value of the constant Co = 2-0r(0. 147088..) = 0.5214.. (see Appendix). In fact, calculating 

the sum of (23) we obtain 0.588.., which is again close to the total value 1/2. This fact is in 
agreement with the small value of the subleading term in the correlator G{x) (see Appendix) 
and the observation [5] that the exact correlator is extremely close to its asymptotic value 
(7) at an arbitrary distance x = 1,...L. The same estimates can be performed for the 
case of the XXZ spin chain which leads to the similar results and in fact can be used to 
obtain the sufficiently accurate estimate for the constant Co [11]. Note, that exactly the same 
relation for the spinless fermion model (Luttinger model) between the constant C(pf, A) in 
the equation (8) and the constant C in the equation (9) takes place: one should simply replace 
the exponent a by a (A). In the context of the possible application for the XXZ model, it would 
be interesting to predict the asymptotic behaviour of the momentum distribution (24) and the 
function (7) from the beheviour of the formfactors (13) for the low-lying excitations (small 
total momentum q). Note also that the calculation of the constant Co in the model of spinless 
fermions (see eq.(8) as well as in the other models (for example, for the XXZ spin chain in 
the magnetic field) remains an open problem. 

Conclusion. 

In conclusion, we presented the new expressions for the formfactors of local operators 
for the XX - quantum spin chain as a Cauchy determinants. Using the functional form of 
the correlator at large distances we proposed the new expression for the constant for the 
asymptotics of the correlator as a Cauchy determinant. Using the scaling arguments, the 
value of the "lowest" formfactor for the XX chain was found. Using the functional form of 
the correlator the momentum distribution for the XXZ spin chain is evaluated. The universal 
character of the connection of the momentum distribution singularity with the asymptotics 
of the correlators and the constants in front of the asymptotics both for the XXZ spin chain 
and the Luttinger model (spinless fermion model) was pointed out. It is possible that the 
formfactor approach can be useful for the calculation of the constant Co in the general case of 
the XXZ - spin model and the other models solvable by the algebraic Bethe ansatz method. 

Appendix. 

For completeness wc present here the formulas for the integrals required for the numerical 
estimate of the constant Ci in the text and present some of the results obtained previously 



15 



for the constant Co- 

First, the integrals one can use are [15]: 



(sin(|/))" ^ V2'2;^V2;V'2 

where B{x;y) = r(x)r(y)/r(x + y), r(l/2) = y^, r(l/4) = 3.625600.., r(3/4) = 1.225417.., 
which is sufficient to estimate the constant Ci = -y/7r(r(l/4) /r(3/4))Co. To obtain the formula 
for the momentum distribution presented in the text one can use the following integral: 



I 

Jo 



1 Tre*^'^/^ 



(sin(y))i-'^ ( 



u+q+l . v-q+1 
2 ' 2 



and take into account the well known relation r(l — x)T[x) = 7r/sin(7rx). The correspon- 
dence between this two formulas can be established with the help of the relation T{2z) = 
(2^^^^7r^^/^)r(2;)r(2; + 1/2). To obtain the asymptotic formula (24) one should use the well 
known expression for the asymptotics of T{z) at large z. 

For completeness let us present here some results obtained previously for the correlator 
G{x) in the thermodynamic limit. Using the expression (18) in the thermodinamic limit we 
calculate the product: 



.ttJ \{2k + l){2k - 1) J 

Considering the logarithm oi Rn after some algebra one can finally obtain the result: 

HR.) = -iln(iV) + i f f (e-« - ^) - (25) 

where the omitted terms arc of order ~ 1/N'^. The above expression coincides with the formula 
proposed in ref.[ll] in the particular case of the XX - chain. The last term in (25) gives the 
following coefficient for the next-to-leading asymptotics for the correlator: 

where the constant Cq is defined in (7). Remarkably, the equation (26) shows that the exact 
correlator G{x) coincides with the leading asymptotic rcsiilt (7) even in the region x ~ 1 
with the sufficiently high accuracy. The value of the constant corresponing to the subleading 
term in eq.(26) was first obtained by McCoy [2] using the method [3] with the help of the 
asympotics of the Barnes G- function [16], defined by G{z -|- 1) = r{z)G{z), G{1) = 1. The 
general expression for the Cauchy determinant is: 



{G{N)yG{l + z)G{l ~ z) 
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In particular, the product Rn can be represented as 

= {Gil/2)) 



(G(7V + i))2 



G'(A^+ l/2)G'(A^ + 3/2) 
Using the asymtotics of the function G{N) at large N, 

G(N) = i - in^ - i(ln2.)iV + (ijV^ - i) In/V - 5iV= + O (ij) . 

where A is the Glaisher constant (see below), one can obtain the result (26) with the constant 
(Co)V2 = 7rV42Vi2ei/M-3. One can evaluate the integral in (25) to get the asymptotic 

HRN) = -\HN)+(^-^ + K'{-i)y 

which is equivalent to the estimate 

R^ = (7V-V4) (2Vi2eV4^-3j ^ A = eVi2-C'(-i) ^ 1.282427 .... 

This result agrees with the result obtained by Wu [3], using the expression of the product 
through the Barnes G- functions [16]. Thus for the constant before the asymptotic (7) the 
known value Co/2-y/7r = 0.147088 ... is obtained. 
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